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ABSTRACT 


We have studied the details of pair annihilation In hot e + -e" plasmas. 
We have calculated the annihilation rate, luminosity and spectrum of 
optically thin plasmas of temperatures above 10 8 K by means of a Monte- 
Carlo simulation. For a given temperature, the spectrum Is peaked at an 
energy equal to 0.511 MeV plus a positive definite quantity of order kT. 

In high temperature sources, such as ganma ray bursts, this blueshift 
can amount to a significant fraction of 0.511 MeV. The annihilation line 
Is also temperature broadened: The width varies as T 1 * for kT much less 
than 0.511 MeV, and as T for kT >>0.51 1 MeV. The widths of the 400 to 
460 keV emission lines observed from several gamma ray bursts, if inter- 
preted as gravitationally redshlfted e + -e” annihilation, set limits on 
the temperatures of the pair annihilation region in burst sources. These 
limits ( T cdxl0 8 K) are considerably lower than the typical kinetic tem- 
peratures (T ?>10 9 K) of the radiating particles. The burst emission, 
therefore, is either nonthermal , or the pair annihilation region is 
spatially distinct from the site of the outburst itself. 

Subject headings: positron-electron annihilation - high temperature 

plasmas - gamma ray bursts. 
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I. INTRODUCTION 

Recent transient gamma ray observations have revealed emission and 
absorption features in many gamma ray bursts. The most commonly observed 
emission line* seen so far < i seven gamma ray bursts (Mazets et al. 1979, 
Teegarden and Cline 1980, Mazets et al. 1«81), and in a longer duration 
transient (Jacobson et al. 1978), falls In the energy range from 400 to 
460 keV. This line could be gravitationally redshlfted e*-e“ annihila- 
tion produced at the surfaces of neutron stars (e.g. Ramaty, Borner and 
Cohen 1973, Mazets et al. 1979, Ramaty et al. 1980). That e + -e“ pairs 

should, In fact, be copiously produced In gamma rc»* burst sources is Im- 
plied by the opacities of the source regions to photon-photon pair 
production (Schmidt 1978, Cavallo and Rees 1978). The exact manner in 
which these pairs would annihilate to produce a relatively narrow line 
from an otherwise optically thick source is not well understood, although, 
as suggested by Ramaty et al. (1980) and Ramaty, Lingenfelter and Bussard 
(1981), synchrotron cooling before annihilation in the strong magnetic 
fields of neutron stars could produce an observable annihilation feature 
in the outer skin layer of the source region. 

Clearly, much information on gamma ray burst sources and neutron 

stars would follow from a complete understanding of the physics of line 
+ 

formation in e -e annihilation. Positron annihilation in relatively cool 

7 

(<10 K) media has been investigated in considerable detail by Crannell 
et al . (1976) and Bussard, Ramaty and Drachman (1979) who applied their 
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results to the observed 0*511 MeV line from solar flares (Chupp et al, 1973) 
and the galactic center {Leventhal, MacCallum and Stang 1978, Jacobson 1981). 

Gamma ray burst sources, however, are likely to have temperatures ccn- 

a 

slderably In excess of 10 k, as Indicated by their observed continuum 

a 

spectra. Whereas below about 10 K positronlum formation plays a central 
role In determining the shape of the 0.511 MeV line (see Ramaty and 
Llngenfelter 1981 for a review), at the temperatures characteristic of 
gamma ray bursts, the bulk of the annihilations proceed directly, 
e + 4-e”'^y+v. Since no intermediate bound state Is produced in such annihi- 
lations, the annihilation feature peaks at an energy which equals the 
electron rest mass, md* =511 keV, plus a positive definite quantity that 
depends on the kinetic energies of the annihilating pairs. This blueshift 
is of the order kT, and thus significant for studies of gamma ray burst 
s’curces, in particular for studies that depend on a precise determination 
uf the gravitational redshift. The fact that the peak of the direct 
e + -e“ annihilation line is temperature dependent and always above 0.511 
MeV, has not been widely appreciated in the published literature. 

In the present paper we calculate in detail the annihilation rate and 
annihilation spectrum of optically thin e + -e~ plasmas of temperatures above 10 a K. 
Our results, while obtained by different kinematical and numerical techniques, 
agree quite well with those of a recent calculation by Zdziarskl (1981). Both 
calculations employ Monte-Carlp simulations for the evaluation of multidimensional 
integrals. We study, in addition, the blueshifting of the peak of the annihi- 
lation spectrum and its broadening as a function of temperature, and we discuss 
the consequences of these effects on the emission lines observed in gamma ray 
bursts. Fhe same calculations could also be applicable to other high energy 
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astrophysical phenomena which Involve radiation from e*-e“ pairs* In 
particular active galactl ■ nuclei and accreting compact x-ray sources. 

While the full understanding of the nature of e + -e" annihilation radiation 
from compact objects will require the study of radiation transfer In a 
relativistic plasma, a necessary Ingredient In such a study Is the anni- 
hilation emisslvlty worked out in the present paper. 

II. HIGH TEMPERATURE e*-e" ANNIHILATION 

Consider an electron-positron pair plasnui In which both kinds of 
particles have equal densities, n + =n - , and isotropic Maxwell-Boltzman 
energy distributions. We shall calculate the annihilation rate and 
luminosity of this plasma, and the energy distribution of the annihilation 
radiation. 

The annihilation rate per unit volume can be written as 

R = n + n_<a(y^,y - , H .)v r (v + ,\.,p)>, (1) 

where cr and v r are the annihilation cross section and relative velocity of 
a positron and electron with Lorentz factors y + and y_, respectively, and 
angle cos"V between their directions of motion. Equation (1) is 
averaged over a uniform distribution for n, and over a Maxwell-Boltzman 
distribution for y + and y_, 

F(v) = (a/K (a))Qy e~ ay , (2) 

2 

a x a 

where p = (y -l) fc 7y, a * me /kT, and 1^ is the modified Bessel function of 


the second kind. 
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To evaluate equation (1) we use the Invariance property of cross 
sections, Y+Y j*v r ■ Invariant (Landau and Llfschltz 1970), and the anni- 
hilation cross section In the center-of-momentum (CM) frame of the two 
particles (Keltler 1954) 


* * ^O 2 

o * — sir it 

4. Y P 



(3) 


a a 

Here r Q * e /me , y* is the posUron or electron Lorentz factor In the CM 

*3 l 

frame and p* * (y -1 )*/y* Thus, 


★ *a 


R = n + n_<a (y )2c0 y /(y + yj>. 


(4) 


where the right hand side is a function of y + , y_ and » only, since 




y* = 2"*[l+y + y.-(y + -l) a (y_ -l)Vf 


(5) 


Equation (5) can be derived from the invariance of the magnitude of the 
total 4-momentum of the positron-electron system. 


We use a Monte-Carlo technique to average equation (4). We choose 
N triads of uniformly distributed random numbers between 0 and 1, x^. , x +i , 

x_.j , and we evaluate , y + ^ and y^ for each triad by solving 
u i = l"2x i and the integrals 

x ±i = /i ±1f (y)<*Y» 


( 6 ) 
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h ft 

where F(y) Is given by equation (2). We then evaluate and a ( y^ ) 
from equations (3) and (5), and the annihilation rate from 


n+n.N 


•’i 

1-1 


y<* 


( 7 ) 


where 


y 1 ■ o*(y*)2c^y’{^Y^V.^). 


(8) 


The results are shown In Figure 1* based on numerical calculations 

with N=5xl0*. At relativistic temperatures R approaches asymtotically 

4xlo' 15 (kT/mc a ) " 8 log (kT/mc 8 ), reflecting the behavior of the 
annihilation cross section at high energies . In the temperature range 

10 a to 10®, R/n+n/nTrJc, because the inverse velocity dependence of the 
nonrelatlvlstic cross section Is cancelled by the relative velocity be- 
tween the particles. At temperatures below 10®K, R/n+n, is taken from the 
calculations of Crannell et al. (1976) and Bussard et al. (1979). Here 
the Coulomb attraction between the particles causes the annihilation rate 
to exceed nr 2 c„ 


The curve L in figure 1 shows the annihilation luminosity. Since the 
energy liberated per annihilation Is mc®(Y + + Y_)» L can be evaluated In 
the same Monte Carlo calculation as R, i.e. 


L s n + n_mc a N 



<v+i + v.,)y i . 


(9) 
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At relativistic temperatures, L also decreases with Increasing T, although 
not as fast; as R. For 10*< r <l(P, L Is approximately constant, because a 
constant amount of energy, 2mc* , Is liberated per annihilation. There Is 
a small positive slope In L(T) around 10 9 K reflecting an Increase In the 
liberated energy without a substantial decrease In annihilation rate. 

The ratio L/2R Is the mean photon energy per annihilation. From the 
numerical calculations we find that for kT«mc 8 , L/2R »• me® ♦ j kT, l.e. 
the mean photon energy equals the electron rest mass energy plus the mean 
particle energy for a nonrelatlvlstlc Maxwell-Boltzman distribution. For 
kT»mc 8 , on the other hand. L/2R 2kT, which Is less than the mean energy 
(3kT) for an ultrarelatlvlstlc Maxwellian. This Is because In the relati- 
vistic regime the annihilation rate decreases with increasing energy causing 
relatively more annihilations of lower energy pairs than of pairs of higher 
energies. In the nonrelatlvlstlc regime, the annihilation rate is energy 
Independent and therefore the mean particle and photon energies are 
essentially equal. 


* 
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To evaluate the energy spectrum of the annihilation radiation, a 
somewhat more complicated approach Is needed since It Is also necessary 
to keep track of the Individual photon energies resulting from each 
annihilation. The positron-electron configuration, specified by 
V+ »Y„ only determines the energy of the two photons in the CM frame, 
twf ■ttog«Y*mc a . The determination of the photon energies In the frame of 
the plasma requires that we also specify the directions of the photons. 

The additional Independent variables are cos’V* and $*, the CM angles, 
respectively, between one of the photons and one of the particles, for 
example the positron, and between the normals to the planes defined by the 
positron and the electron and the positron and the gamma ray. The distribu- 
tion of Is uniform, while that of ^* Is determined by the differential 
’iross section in the CM frame (Heltler 1954) 


da* 


nrfl 1 T )(2-u 

-r 


asfl 



The photon energies and fiufe In the frame of the plasma can be 
evaluated as follows: Let ? c be the velocity of the CM frame, and 
cos” 1 ^*^ c and cos” 1 n* +c be the CM angles between one of the photons and 


wm 
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3 C and the positron and ? c , In terms of the total energy and momentum 
of the positron and electron, E++E_ and ? c and the corresponding 
Lorentz factor y c » can be written as 

t e ■ C(P + +r. )/(E + +E_). _ (II) 


and 


Y c * (Y+ + Y. )/(2 y*) • 02) 

where we use equation to derive equation (12) f'om equation (11). 
Then and are given by 

fw»j * me 3 y* Y c ( 1 =*=0 yc ^ * 03) 


where 


n* c * 4 C »M1<)^0-^)W 


04) 


and 


Vv c i* 

<’i*-1)V a -0 ,i 05) 


Thus, given y + * Y.» u>* and $*, we can evaluate the photon energies In 

the plasma frame by using equations (5), (12), (13), (14) and (15). As a 
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consistency check, we verify energy conservation: from equations (12) and 
(13) It follows that tkoj ♦Ktfe»mc 8 {Y 4 *Y.)» 


In the numerical evaluation of the energy spectra we follow a similar 
procedure as In the evaluation of the annihilation rates. Since we now 
have 5 Independent variables, we choose two additional uniformly distributed 
random numbers between 0 and 1, x^* ?,nd x^, evaluate v + ^, y.j and ^ as In 
the calculation of R, calculate yJ B l-2x # and $?*2ttx.*,‘ and then evaluate 

» yj i 'fj 

the quantity 


AR(E) « n+njf 1 ga*(vj)]“ l 2da*/dy*(Y**y*)y^ * 


(16) 


where is given by equation (8), For a given photon energy, E, the 
summation Is over all values of 1 for which either ^ or fagj, determined 
from equations (5), (12), (13), (14) and (15), are in the range E to E+oE. 
In the numerical calculations, we choose AE sufficiently small to obtain 
an essentially continuous spectrum. 


Numerical results for N*5xl0 6 are shown In Figures 2 and 3, where 
(2/n + nJ dR/dE is plotted as a function of photon energy for several 
plasma temperatures from 10 a to 10* 1 K. As expected, the peak energy of 
the annihilation line, E m , and its full width at half maximum, FWHM, 
Increase with T. At temperatures much larger than mc a , the annihilation 
spectrum bears no relationship to tne characteristic energy 0.511 MeV; 
instead it is determined by the convolution of tne energy spectrum of the 
e* -e’ pairs and the annihilation cross section. 


&»»55e3r?fei 




WMSswsw 



10 


By using the numerical results as shown <n Figures 2 and 3, and by 
evaluating annihilation spectra also for a number of other fcsnperatures, we 
have deduced the dependences of E^-mc 8 and the FWHM on temperature* Me show 
these in Figure 4, For kTc^c 3 , we can approximate the width by the non- 
relativistic formula FWHM " v lleV(T(K)) v (e.g., Crannell et al. 1976) and the 
blueshift by E^mc 3 (3/4)kT» For kTSMitc 8 both the width and the shift vary 
linearly with T, in particular £ w °- l,2kT. By comparing these results with 
the mean photon energy discussed above, we see that the peak of the annlhlla- 
t*un spectrum is always at a lower energy than its mean* a result caused by 
the asymmetry of the photon distributions, as can be seen, for example, 
in Figure 2. 

We have repeated the calculations of the present Section using a 
Fermi-nirac distribution with zero chemical potential, F(yHpY s /(e a *+1 ) , 
instead of the Maxwellian given by Equation (2). This distribution 
describes e + -e" pairs in equilibrium with black-body radiation (e.g. 

Landau and Lifshits 1958). Without considering in detail which distribu- 
tion might be more appropriate for gamma-ray burst sources, we note here 
that the two distributions differ appreciably only for kT £ me 3 . Indeed, 
the annihilation rate turns out to be smaller for the Ferroi-Dirac case 
than for the Maxwell -Boltzman case by only about 18 at 3xlO a K, 8 % at 
lff J K and 15" at 10 n K. The blueshifts are also not very different for 
these temperatures. At lO^K, for example, E rn is larger for the 
Fermi -Dirac distribution by only about 10" than the E m for a Maxwellian 
at the same temperature. As we shall see in the next Section, the 
temperatures of the annihilation regions of gamma-ray bursts must be 
quite low 3xlO°K) to produce lines of widths that do not exceed the 
observations* At these temperatures the differences between the two 
distributions are essentially negligible* 
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III, DISCUSSION AND CONCLUSIONS 

The results of our calculations could be applied In a straightforward 
manner to observations of annihilation radiation from optically thin sites. 
Such sites might be present, for example. In the vicinity of active 
galaxies, although observations exist only from the GalacMc Center where 
the annihilation site is cooler than 10 6 K (Bussard, Ramaty and Drachman 
1979). 

Gamma ray burst sources are likely to be optically thick. Neverthe- 
less, we can apply our results to their last optical depths, which must 
be the regions in which the e + -e” pairs responsible for the observed 
emission lines actually annihilate. An upper limit on the temperature of 
these annihilation sites can be obtained from the observational upper 
limits on the line widths. For the 400 to 460 keV lines, which we assume 
to be gravitationally redshifted e + -e“ annihilation, the full widths at 
half maximum are less than 200 keV (Teegarden and Cline 1980, Mazets 

et al. 1981). Then from Figure 4, the temperature of the annihilation 
region is less than 3x1 0 8 and hence E m -mc 2 <20 keV. This shift, about 4% 


In the line energy, Is significant, since the surface redshifts of 
neutron stars range from about 5 % to 50% depending on their mass and 

equation of state (Borner and Cohen 1973). 

We should note, however, that a blueshlft that does not exceed 
20keV is actually smaller than the observational uncertainty In the peak 
energies of the lines (Teegarden and Cline 1980). Also, In the deter- 
mination of a neutron star parameter, such as mass, the uncertainty 
caused by a blueshift ranging between 0 and 20keV Is of the same order 
or less than the uncertainty Introduced by the poor knowledge of the 
neutron star equation state (see calculations of Borner and Cohen 1973). 
Nevertheless, future spectroscopic observations with higher sensitivity 
could determine the peak energy of the line and Its width much more 
accurately, and these could provide valuable information on neutron star 
equations of state (see for example Brecher 1977). 

Tho upper limit of 3xlO a K on the temperature of the annihilation 
region should be compared with the energies of the particles that produce 
the observed continuum emission of the bursts. Since all gamma ray bursts 
with observed 400 to 460 keV emission lines have continuum spectra which 
extend to at least several hundred keV (Mazets ef. al. 1981, Teegarden and 
Cline 1980), the radiating particles must have energies of at least these 
values, i„ec higher by an order of magnitude or more than the kT of the 
annihilation region. Furthermore, conditions imposed by the pair produc- 
tion threshold, also require particles or photons of energies at least 
comparable to rnc a to produce sufficient pairs . Thus, either the positrons 
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annihilate In a spatially distinct region from that responsible for the 
continuum and the pair production, or the annihilation region Is non- 
the.mal . 

A spectral feature with peak energy around 430 keV has also been 
observed from the March 5, 1979 transient (Mazets et al . 1979). The 
temperature-induced blueshlft, discussed In the present paper, was not 
taken Into account In the neutron star vibrational model of this transient 
proposed by Ramaty et al. (1980). These authors relate the duration of 
the impulsive spike of the transient to the gravitational redshlft by 
arguing that the duration could equal the damping time of the vibrations. 

We find here that with an upper limit of 20 keV on the magnitude of the 
blueshift, and considerable uncertainties in both the energy of the line 
peak and equation of state, these arguments are not affected by the 
neglect of the blueshift. Regarding the difference between the radiation 
and annihilation temperatures, the neutron star vibrational model is 
expected to produce nonthermal emission, as discussed by Ramaty et al. (1980) 
and Ramaty, Llngenfelter and Bussard (1981). 

As already mentioned in the Introduction, the blueshifting of the 

e + -e" annihilation line in hot plasmas is the consequence of tha fact 

that direct e + -e" annihilation proceeds without the formation of an 

intermediate bound state. On the other hand, a gamma ray line resulting 

from the deexcitation of a bound nuclear state, for example the line from 
12 

excited C, peaks at 4.44 MeV, independent of the temperature of the 
plasma. The production of this line in accreting matter around neutron 
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stars and black holes has been considered by Higdon and Llngenfelter 
(1977). However, In e + -e* annihilation, the only available bound state 
Is posl tronium, and at the high temperatures of interest here, the rate of 
posltronium formation by radiative combination, e + +e" + Ps+y, Is negligible 
in comparison with direct annihilation e + +e’-*-y+Y. 

Another gamna ray line of considerable astrophysical Interest which 
is formed without an intermediate bound state is the 2.221 MeV line from 
neutron capture on hydrogen, n+ 1 H-* 2 H+y (e.g. Ramaty, Kozlovsky and 
Lingenfelter 1975). In the 20-minute gamma ray transient of Jacobson 
et al . (1978), one of the observed lines has been identified with neutron 
capture on 1 H. As for the 0.511 MeV line, for this line as well, the 
temperature-induced blueshift is of the order kT. We defer a more de- 
tailed evaluation of this shift and its astrophysical consequence to 
another paper. 

In conclusion, we have pointed out a hitherto ignored effect iii 
astrophysical applications, the blueshifting of the e -e" annihilation 
line due to the temperature of the annihilation region. We have evaluated 
in detail the annihilation rate and annihilation spectrum of a hot e -e 
plasma. By applying our results to gamma ray bursts, we find that the 
temperature of the annihilation region in burst sources is less than about 
3x1 0 8 K, that the temperature-induced blueshift of the 0.511 lieV line is 
less than 20 keV, and that either the annihilation region is nonthermal 
or spatially distinct from the burst site itself. 
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FIGURE CAPTIONS 


1. Annihilation rate, R, and annihilation luminosity, L, of an e + -e' 
plasma as functions of temperature. 

2. Annihilation spectra for various temperatures. 

3. Annihilation spectra for various temperatures. 

4. Temperature-Induced blueshlft, E m -mc 2 » and full width at half 
maximum, FWHM, of the e + -e“ annihilation feature as function of 
temperature. 
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